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1.‘  Introduction.  The  classical  noisy  duel,  say  of  1  bullet  vs  2, 
with  equal  accuracy  and  a  continuum  of  opportunities  for  firing,  does  not 
have  a  saddle  point.  Remedies  have  included  1)  "continuizing  these  ammuni¬ 
tion  stores  ([12],  [ 14 ] ) ,  2)  settling  for  only  e-good  solutions  ([2],  [11]), 
and  3)  placing  certain  extreme  behavioral  demands  ([4])  on  the  weaker  duelist, 
calling  for  anticipating  the  opponent's  current  action,  i.e.,  for  instanta¬ 
neously  adapting  his  own  firing  plans  to  that  action.  A  fourth  possible 
remedy,  introduced  in  [15]  and  in  a  sense  dual  to  the  first,  is  to  allow 
only  a  finite  set  of  firing  opportunities,  thus  creating  approximating 
finite  simultaneous  ga.nes  ([17]),  which,  as  finite  games  of  perfect  recall 
([13]),  admit  randomized  behavioral  solutions,  i.e.,  solutions  in  the  form 
of  firing  time  distributions.  The  asymptotic  versions  of  these  distribu¬ 
tions  model  the  anticipation  required  of  the  weaker  duelist  in  [4]. 

A  not  unimportant  element  of  zero-sum  two-person  game  theory  is  the 
quantification  (typically  in  stochastic  fashion)  of  compelling  not  other¬ 
wise  quantified  features  of  strategy;  the  quantification  [19]  of  bluffing 
in  poker  for  example,  by  behavioral  randomization  calling  for  "irrational" 
play  now  and  then;  or  the  quantification,  by  an  atom  of  firing  probability 
reserved  for  the  bitter  end,  of  the  Kamiki^e  ommitment  of  a  value  -  less 
single  -  salvo  fighter  in  silent  duel'  wit.'  vtinuously  firing  opponent. 

Still  another  example,  then,  is  our  modeling  of  the  weaker  duelist's 
instantaneous  anticipation  by  the  fact  that  his  asymptotic  first-firing 
distribution  has  infinite  expectation  and  straddles  that  of  his  opponent. 

It  seems  fair  to  say  that  this  modeling  will  be  of  interest  to  the  degree 
that  its  details  will  not  be  predictable  intuitively  a  priori,  and  in¬ 
deed  they  do  not  seem  to  be.  It  does  seem  impossible  to  foretell,  start¬ 
ing  only  from  knowing  that  the  players'  good  strategies  somehow  are  to 
be  asymptotically  stochastically  portrayed,  whether  the  same  normalization 
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will  work  for  both  players,  where  the  two  asymptotic  firing  time  distri¬ 
butions  will  be  situated  with  respect  to  one  another,  or  whether  I's  need 
to  anticipate  wilj.  be  reflected  in  incomplete  convergence  in  law  of  his 
normalized  firing  time,  beyond  the  infinite  expectation  of  its  asymptotic 
distribution.  What  does  turn  out  is  given  below  in  (1) .  Whether  portions 
of  (1)  are  extrapolatable  from  T(l,2)  to  games  of  timing  generally,  we 
do  not  as  yet  know,  although  we  do  know  that  they  are  so1  extrapolatable 
to  r(m,n),  as  outlined  in  Section  6. 

Consider  then  two  equally  accurate  duelists,  players  I  and  ,11,  possess¬ 
ing  respectively  1  and  2  bullets,  facing  a  finite  grid,  with  uniform  mesh 
e,  of  joint  firing  opportunities,  with  the  points  of  the  grid  calibrated 
from  0  to  1  according  to  bullet  lethality  p.  Dynamic  programming  yields 
an  optimal  behavioral  solution  for  each  duelist,  which,  among  other  things, 
specifies  optimal^probabilities,  greater  than  zero  and  less  than  one,  of 
first  firing  (strictly  speaking  just  "firing"  in  the  case  of  I)  at  certain 
contiguous  grid  points  near  p  =  1/3.  These  probabilities  are  conditional 
on  neither  duelist  having  yet  fired,  and  may,  for  each  duelist,  be  inter¬ 
preted  as  determining  an  optimal  unconditional  distribution  of  time  of 
(first)  fire  near  p  ■  1/3,  say  <J>1  for  player  I  and  4>  for 
player  II.  Players  I  and  II  are  to  Independently  select  the  time  of  their 

first  shot,  respectively  according  to  <!>  and  *f>  ,  and  to  fire  at  the 

1 1 C  Z  i  c 

time  selected  unless  the  opponent  has  fired  before  that  time. 

Though  both  41-  and  ‘1>  „  converge  to  degeneracy  at  p  =  1/3 

1  >  C  L  ,  C 

with  increasing  numbers  of  firing  opportunities,  the  microscopy  of 

location-scale  norming  (the  same  norming  for  4’  and  <!'„  )  reveals 

I  I ,  £  — *  £ 

two  asymptotic  distributions  on  [-!,+»),  say  <J>  and  0^,  with  densities 
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As  indicated  in  Section  5,  relation  (la)  may  be  derived  with  the  aid  of 
certain  relatively  weak  hypotheses  clearly  indicated  by  extensive  computer 
runs;  relations  (lb),  (lc)  and  (Id)  follow  from  certain  strengthened  ver¬ 
sions  of  these  hypotheses,  allowing  the  expression  (34)  of  0^  and  (38) 
of  »i>2  *n  terms  of  certain  Hermite  equation  solutions. 

We  submit  that  (la)  and  (lc)  model  I's  need  to  anticipate  by  I's 
asymptotic  firing  distribution  "straddling"  that  of  II, ''and  that  this 
need  to  anticipate  is  also  modeled  by  the  upper-tail  mass  of  I's  asymptotic 
firing  distribution,  mass  enough  to  cause  (Id)  (though  not  enough  to 
prevent  (lb)). 

The  next  4  sections  detail  the  origins  of  the  findings  (1)  concern¬ 
ing  ^  and  $2*  together  with  the  derivation  of  a  certain  associated 
differential  equation  whose  solution  is  discussed  in  [9].  Section  6 
indicates  the  extension  i>f  che  results  to  f’(m,n). 

2.  The  dynamic  programs.  As  pointed  out  in  the  introduction, 
restricting  firing  opportunities  to  a  finite  grid  converts  the  noisy 
duel  into  a  finite  simultaneous  game  (cf.  [17]).  The  successive  infor¬ 
mation  sets  of  this  simultaneous  game  corresponding  to  no  one  having 
yet  fired  tell  the  story:  They  are  the  only  non-trivial  information 

v 

sets  for  the  duel  T(l,2)  treated  here.  Since  each  player  possesses 
precisely  one  such  information  set  per  grid  point,  we  shall  speak  below 
of  a  player's  behavior  "at  a  grid  point  p",  when  we  in  fact  mean  to 
speak  of  his  behavior  at  the  information  set  corresponding  to  that  grid 
point. 

Assuming  a  uniform  grid  of  mesh't,  at  each  of  the  successive  grid 
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points  p  **  1,  1-e,  l-2e,...,0  proceeding  backwards  from  unity  the 
players  face  a  2x2  game  with  matrix 


(2a) 


(2b) 


* 

F 

F 

'  "  L  ' 

0 

1 

? 

-1 

0 

F 

F 

F 

-(1-P)2 

2p-l 

?' 

-P 

u[p+e] 

when  p  =  1,  and 


when  p  <  1, 


where  u[p+e]  is  the  value  of  the  game  at  the  grid  point  p+e  (see 
[18]).  Since  2p-l  is  never  the  minimum  in  its  row  for  p  >  0,  the 
game  has^ a  saddle  point  at  (F»F)  as  long  as  p  p*,  where  p*  »  0.382 
is  the  solution  of  -  /p*. 

At  a  grid  point  p  <  p*,  one  of  three  things  can  happen.  There 
may  be  a  saddle  point  at  (F,F),  i.c.,  -p  u[p+e]  _>  2p-l ,  which  forces 
P  <  1/3. 

Also,  there  are  grids  of  arbitrarily  small  mesh  for  which  some 
grid  point  p  <  p*  produces  a  saddle  point  at  (F,F),  For  example, 
if  <S  >  0  and  N  is  a  positive  integer  such  that 

2N-1  >  max{2(l-p*)/6,  (1-p*) /(l-2p*) } , 


then  the  grid  of  mesh  c  °  2(l-p*)/(2N-l) ,  which  is  less  than  5, 
produces  a  (F,F)  saddle  point  when  p  -  1-Ne  ■»  p*~e/2.  It  is  shown 
following  (A)  that,  for  given  e,  a  ^F,F)  saddle  point  can  occur  at 
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most  ond  grid  point  p. 

The  third  possibility  for  a  grid  point  p  <  p*  is  that  the 
corresponding  game  has  no  pure  saddle  point.  As  is  shown  in  [18],  the 
value  u [ p+c ]  determines  the  mixed  value  u[p]  and  the  optimal 
conditional  firing  probabilities  x^[p]  for  I  and  x^lp]  for  II: 

(3a)  u[p]  =  ~{u[p+e] (1-p) 2  +  p(l-2p)}/{u[p+e]  +  p(l-p)), 

(3b)  x^p]  ®4u[p+e]  +  p}/{u[p+c]  +  p(l-p)}, 

(3c)  x2[p]  =  {u[p+c]  +  l-2p)/{u[p+e]  +  p(l-p)}. 

By  inspection  of  the  matrix  (2b)  of  the  game,  we  have  in  turn  that 
-  (1-p)2  <  -p,  i.e.,  p2  -  3p  +  1  >  0,  that  u[p+e]  -p,  and  that 
u[p+c]  <  2p~l.  The  last  inequality  implies  that  the  denominator  in 
each  o£_the  formulas  (3)  is  less  than  -p2  +  3p-l,  which  is  negative. 
Evidently,  then,  both,  x^fp]  and  x2[p)  are  positive.  Moreover,  it 
is  straightforward  to  show  that  formulas  (3)  imply  that 

(4a)  u[p]  +  p  =  -  (p2-3p+I)  x^p], 

(4b)  u[p]  +  l-2p  =  -  p2  x2fp], 

(4c)  u [ p+e ]  -  u[p]  =  {u[p+e]  +  p(l-p)}  x x [ p ]  x2[p), 

A. 

all  of  which  are  negative. 

In  particular,  u[p)  <  -p  <  -(p-e),  so  the  game  at  the  next  grid 
point  p-c  has  no  saddle  point  at  (F,F).  This  implies  that  if,  in 
the  backwards  induction  process,  a  (F,F)  saddle  point  occurs  for  the 
first  time  at  a  grid  point  p+e,  where  e  <  (p*-l/3)/2,  then  necessarily 
p*-e  <  p+e  _<  p*.  Also,  if  this  happelvs,  then  the  game  for  the  next  grid 
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point  p  has  the  matrix  (2b)  with  u[p+e]  =  -(p+c),  so  it  does  not  have 
(F,F)  as  a  saddle  point.  This  discussion  shows  that  the  (F,F)  saddle 
point  phenomenon  pan  occur  at  most  once  for  a  given  mesh  e,  either  at 
p*  (if  that  happens  to  be  a  grid  point)  or  at  the  first  grid  point  less 
than  p*. 

Hence,  we  can  conclude  that,  for  e  small,  there  is  a  largest  grid 
point  p(e)  at  which  the  game  has  no  pure  saddle  pointy  and  that 
p*  -  2 c  <  p(e)  <^p*.  Starting  at  p(c)  and  continuing  backwards,  we 
move  through  a  succession  of  grid  points  p  at  which  mixed  optimal 
strategies  obtain.  Since  (4c)  implies  that  u[p+e]  <  u[p],  the  value  of 
the  game  increases  as  we  proceed  backwards.  Moreover,  the  game  at 
P  =  p(e)  +  £  has  a  pure  saddle  point  at  either  (F,F)  or  (F,F).  So, 
the  value  u[p(e)  +  e]  has  one  of  the  forms  -(1-p)2  or  -p,  and  is 
therefore  greater^ than  M.  It  follows  that  we  must  eventually  arrive 
at  a  smallest  grid  point  p(e)  for  which  u[p(e)  +  e]  <  2p(c)  -  1,  i.e., 
for  which  random  strategies  occur. 

Finally,  the  matrix  (2b)  for  the  game  at  p  =  p(e)  -  e  has  u[p(e)] 
in  the  (F,F)  entry,  so  it  has  a  saddle  point  there.  Proceeding  induc¬ 
tively  backwards,  then,  the  same  situation  holds  for  0  <_  p  <  p(c). 

The  conditional  firing  probabilities  x^[p(e) ] , . . . ,x^[p(c) ]  deter¬ 
mine  the  unconditional  firing  time  distribution  <J>  of  the  introduction. 

An  analogous  interpretation  holds  as  well  for  the  sequence  x^IiKe)  ] , . . .  .x^tp^)  ] 
and  the  corresponding  distribution  4>  .  The  next  sections  discuss  the 

scaling  of  •  and  ,  and  consequent  derivation  of  their  asymptotic 

1  ,  C  L  ,  C 

counterparts  and  ( 
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3.  ^Asymptotic  scaled  values  and  probabilities  of  fire.  For  a 

finite  grid  (1,  1-e,  l-2c, . . . ,0) ,  let  p(e)  and  p(e),  respectively, 

be  the  smallest  and  largest  grid  points  from  Section  2  at  which  the 

corresponding  2x2  game  has  no  pure  saddle  point.  Also,  let  8(c)  =  l/3-p(e), 

let  0(c)  ■  e/8(c)  and,  for  t  _>  -1,  let  p(t,e)  be  the  largest  grid 

point  that  does  not  exceed  1/3  +  tB(e).  This  allows  writing 

*  '• 

(5a)  ^  p(-l,c)  =  p(e) 

and,  for  t  >  -1, 

(5b)  p(t,c)  =  1/3  +  t3(e)  +  f(t,c), 

where  -  e  <  f(t,e)  £  0.  We  introduce  the  step-functions  u(t,c)  =  u[p(t,c)), 

x1(t,e)^  =  x^[p(t,c) ] ,  x2(t,c)  =  x2[p(t,e)),  and  note  that  the  definitions 
of  p(t,c)  and  O^c)  imply  that  u(t  +  0(c), c)  *  u[p(t,e)  +  C]  and 
similarly  for  the  other  two  functions  evaluated  at  (t  +  0(e), c) . 

Extensive  computer  implementations  in  [18]  of  the  dynamic  programs 
of  Section  2  point  to  the  following  hypotheses  concerning  B(c)  and 
u(t,e) : 

HYPOTHESIS  A.  There  is  a  positive  number  C  such  that 

V 

(6a)  8(e)  »  c/F  +  o(/e). 

HYPOTHESIS  B.  There  is  a  function  v(t) ,  differentiable  for  t  >  -1 
and  right-differentiable  at  t  =  -1,  such  that,  if  t  2.  “1*  then 

u(t,c)  -  -  1/3  -  v(t)0(e)  +  g(t,e), 

I  . 


(6b) 
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where  for  each  t,g(t,e)  =  o(/e)  and  g(t+U(e),c)  -  g(t,e)  =  o(c). 

A  simple  application  of  the  definition  of  the  derivative  v(t)  results  in 

(7)  u.(t+G(e),e)  -  u(t,e)  =  -  v'(t)0?(e)  +  h(t,e). 


where  for  each  t,  h(t,e)  =  o(e).  We  proceed  to  use  results  from  Section 
2  to  derive  a  differential  equation  that  v(t)  must  satisfy  under  Hypoth¬ 
eses  A  and  B.  ■' 

P- 

First,  for  t  _>  -1,  we  can  choose  e  sufficiently  small  so  that  the 
2x2  game  for  the  grid  point  p(t,e)  has  no  pure  saddle  point,  i.e.,  so 
that  p(e)  £  p(t,e)  £  p(e).  This  holds  by  definition  when  t  =  -1.  For 
t  >  -1,  the  inequality  p(c)  £  1/3  +  t3(e)  is  equivalent  to  0  £  (1+t) 3(e) //E". 
But  this  must  hold  for  e  small,  since  Hypothesis  A  implies  that  as 
c  -*  0+,  (1+t) 6(e) //£*-»■  (l+t)C,  which  is  positive.  Also,  1/3  +  t3(e)  -*•  1/3 
while  the  inequalities  p'*-2e  <  p(e)  <  p*  force  p(t)  -*•  p*  as  c  -*  0+, 
so  1/3  +  tG(e)  <  p(e)  for  small  e. 

Hence,  we  can  apply  (4c)  to  (7),  obtaining 


(8)  -  V'(t)  +  =  [u(t+0(c),E)  -  u(t.e)l/02(£) 

02(e) 

x1(t,r,)  x2(t,e) 

=  {u(t+0(e),e>  +  p(t,e)(l-p(t,c))}-~^ - — 


The  definition  of  9(c),  together  with  Hypothesis  A  and  the  fact  that 
h(t,e)  «*  o(c),  imply  that  h(t,e)/02(c)  0  as  e  -*■  0+.  And  relations 

(5)  and  (6)  imply  that  the  expression  in  braces  on  the  right  side  of  (8) 
approaches  -  1/9  as  e  ■*  0+. 

Relations  (4a)  and  (4b)  can  be  used  in  conjunction  with  (6b)  to 


calculate 


(9) 


x  (t)  -  11m  x  (t,e)/0(e)  for  1  =  1,2. 
1  .  e+0+  1 


For,  if  we  use  (4a)  we  can  write 


(10)  + 

tV  ’  ' 

If  we  use  relations  (5)  together  with  Hypotheses  A  and  B,  we  obtain,  as 


c  -*•  0+, 


(ID 


xx(t)  =  9[v(t)  -  C2t], 


Similarly,  if  we  start  from  (4b),  we  get 


(12) 


x2(t,e)^ 

6(c) 


Mt. 0+1/3  2(l/3-p(t,Q) 

\  0(e)  +  0(e) 


*  P  (t,e) , 


which,  as  €  -*  0  ,  yields 


(13)  x2(t)  =  9[v(t)  +  2C2t). 

Returning  to  relation  (8)  with  these  results  and  letting  e  -*•  0+,  we  find 
that  for,  t  _>  -1,  v(t)  must  be  a  solution  of  the  Riccati  equation 


(14)  v'(t)  =  1/9  x2(t)  x2(t)  =  9[v(t)-C2t][v(t) +2C 2 1 ] . 

We  can  obtain  an  initial  condition  for  v(t)  at  t  =  -1  by 
examining  (12)  with  t  =  -1  or,  what  amounts  to  the  same  thing,  (4b) 
with  p  »  p(e).  Noting  that  0(e)  >  0  for  e  small  (since  1/C 

as  e  -*•  0+),  we  have  from  (4b)  that 


11 


i 


(15)  ‘  0  >  -p2(c)x2[p(e)]/0(c)  »  (u[p(e)]+l-2p(e))/6(c). 

On  the  other  hand,  the  definition  of  p(c)  and  consideration  of  the 
2x2  matrix  (21)  for  p  =  p(t:)-e  shows  that  u[p(c)]  2  2(p(c)-c)  -  1. 
Hence,  the  right  member  of  (15)  is  not  less  than  -2r./0(e)  =  -2(i(c) , 
which,  approaches  0  as  t:  -*•  0  .  Since  p2(t')  -*•  1/9  as  r.  -*•  0  ,  we 
conclude  that 

(16)  f  x  (-1)  =  lim  x  (-l,e)/6(c)  =0. 

c-*-0+ 

By  (13),  this  can  be  written  as 

(17)  v(-l)  =  2C2 . 

Furthermore,  the  limit  relations  (9)  together  with  the  negativity 
of  the  quantities^!/*  (4a>  and  (4b)  imply  that  x^(t)  _>  0  and  x2(t)  _>  0 
t  _>  -1,  Hence,  (14)  implies  v’(t)  2  0.  Applying  this  to  (13)  gives 
x^(t)  2.  l8c2  >  0.  Hence,  by  (16), 

(18)  x2(t)  ■>  0  for  t  >  -1.  __ 

To  show  that  x^(t)  is  strictly  positive,  suppose  x^(t^)  »•  0  for 
some  tg.  Then  (11)  and  (17)  rule  out  t^  =  -1.  But  if  this  were  to 
occur  at  some  t^  >  -1,  then  the  fact  that  x^(t)  2.  0  for  t  2 
would  imply  that  x^(t)  would  have  a  global  minimum  at  an  interior 
point  of  its  domain.  So,  ( ^q)  “  0  which,  by  (11),  says  v'(tg)  =  C2. 
But  (14)  with  x^(tg)  -  0  says  v'(tg)  ■  0,  a  contradiction.  Therefore, 

x^'(t)  >  0  fok  t  2  *!• 


(19) 
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This,  coupled  with  (11)  and  (12),  gives 


X2<t)  >  27C2t  for  t  _>  -1, 


which  is  a  stronger  statement  than  (18)  for  t  >  0. 

J.  C.  Kegley  has  proved  in  [9]  that,  given  any  C  >  0,  there  is  a 
unique  solution  of  (14)  that  satisfies  (19)  for  large  values  of  t,  and 
that  the  solution  exists  for  ail  t.  Two  forms  of  the  sblution  are  given 
in  terms  of  the  confluent  hypergeometric  functions.  One  of  the  forms, 
which  will  be  used  in  Section  5,  is 


v(t)  .  C2t 

v(t)  L  1  9  <Ks)  ’ 


where  6  =  C/27/2,  s  =  6t,  and  iKs)  is  the  solution  of  the  Hermite 
equat  ion"  - ... 


(22a) 


f'(s)  -  2s  ^’(s)  -  2/3  i^(s)  =  0 


that  satisfies  the  initial  conditions 


ip(0)  =  1,  (0)  =  -2r(2/3)/r(l/6) 


The  second  form  of  the  solution  is  used  to  prove  that  condition  (17) 

V 

then  determines  C  uniquely,  namely,  C  =  -  z^l/JT,  where  z^  is  the 
zero  of  the  Weber  parabolic  cylinder  function  U2/3^'  A  nun,erical 
approximation  of  z^  furnished  by  R.  J.  Lambert  gives 


(22c) 


C  -  0.09CG37  +  3x10  . 


This  is  in  good  agreement' with  the  vafnc  .09067  of  P(r)//c"  obtained 


for  e  ^  1.5x10  the  smailest  of  the  e's  underlying  the  computer 
runs  leading  to  Hypotheses  A  and  B,  and  enhances  the  credibility  of  these 
hypotheses. 

4.  Asymptotic  distributions  of  first  fire.  Consider  now,  for 
t  _>  -1  and  e  small  and  positive,  the  probability  <1>  (l/3+2(e)*t)  H 

!-<!>.  (l/3+S(£)*t)  that  player  1  fires  at  or  after  the  grid  point 

-i 

p(t,e)  +  c  given,  as  always,  that  player  II  has  not  already  fired.  If 
we  denote  by  I(t,e)  the  collection  of  all  grid  points  p  with 
p(e)  £  P  £  p(t,e),  where  e  is  small  enough  so  that  p(e)  £  p(t,e)  £  p(e), 
then  we  have,  by  independence  of  random  moves  under  a  behavioral  strategy, 

(23)  7.  (l/3+P(£) *t)  =  H  (l-x.[p]). 

pCl(t.G)  1 

By  analyzing  the  natural  logarithm  of  4>  (1  /  3+6(0  •  t)  ,  we  will  show 

i  f  c 

that,  as  c  -*■  0+,  7  (l/3+6(c)-t)  tends  to  the  limit 

1  »  t 

(24)  ^(t)  =  exp(-  xx(T)dT)  . 

Thinking  of  t  2  *1  as  fixed,  recall  now  that,  as  indicated  above 
following  (5b),  x^(T,e)/0(e)  is  a  step-function  for  -  1  _<  T  t,  with 
uniform  step-size  0(e).  Hence,  we  may  write 

(25)  Z  x.lp]  =*  /C(E)(*  1(T,e)/0(e)]dt, 
pcl(t.e) 

where  t(e)  t  <  t(e)  +  0(e)  by  definition  of  p(t,e). 

In  order  to  show  that,  as  £  ■*  0+,  x^(x,L)/0(e)  tends  to  x^(t) 

I  . 

uniformly  for  -1  £  T  t,  we  invoke  Hypothesis  B.  This  hypothesis  says 


I 
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that  the' step-functions  (u(t,e)  +  l/3]/0(e)  converge  pointwise  to 
which  is  assumed  to  be  differentiable  and,  hence,  continuous. 

But  (4c)  and  the  fact  that  Q(c)  >  0  shows  that  for  each  e,  the 
corresponding  step-function  is  a  non- increasing  function  of  r.  By 
a  well-known  theorem  (see  for  example  page  90  of  [3]),  it  follows  that 
(u(T,e)  +  1/3 ]/0(e)  converges  uniformly  to  -v(r)  on  the  compact 
interval  f-l,t].  '* 

Inspection  of  relation  (10)  and  the  definition  of  p(t,c)  now 
makes  it  clear  that  the  step-functions  x^(t,g)/0(c)  converge  uniformly 
to  x^(t)  on  [-l,t],  which,  since  t(e)  -*•  t  as  c  ■*  0+,  allows 
concluding  from  (25)  that 

(26) .  T.  x  ( p ]  -*•  /C.x  (i)dr  as  c  ■*  0+  . 

pcl(t.c)  1 

•  -  A 

We  also  note  that,  since  G(r.)  ->-0  as  C  -*•  0+,  the  uniform  convergence 
of  (l ,e)/0(e)  to  x^(t)  for  -1  <  t  <  t  forces 

(27)  sup  x.[p]  -*■  0  as  £  •*■  0+. 
pcl(t,e)  1 

Now  we  examine 

(28)  log  4>  (l/3+0(e)’t)  +  E  x  Ip]  -  l  {log(l-x  [p])+x  [p]}. 

i,e  pel(t.c)  L  peKt.c)  1  1 

By  Taylor's  theorem,  we  have  log(l-x^[p])  +  x^(p]  *  r(p),  where 

J  r  (p)  i  (x1[p])2/,2(l-x1(p])2.  If  we  choose  c  small  enough  so  that, 

by  (27),  sup  x.[p]  £  1-1//2,  we  find  that  |r(p)|  £  (x^p))2  for 
pcl(t.f)  1  I,  ~ 

prl(t,r).  Therefore  the  right  member  of  (28)  has  absolute  value  not 


exceeding 
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1  E  (x.[p])2  £  sup  x  (p]  •  E  x[p], 

pel(t.c)  pd(t,e)  pel(t,c) 

which  approaches  0  as  c  -*■  0  by  (26)  and  (27).  By  (26),  then, 

log  4^  c(l/3+8(e)  •  t)  -*■  -/^x^dt  as  c  ■*  0+, 

so  <J  (l/3+3(c) *t)  converges  to  <1>  (t)  as  defined  by  (24). 

*  #  c  x.  ,t 

Hence,  4*,  (l/3+3(e)*t)  =  1  -  4>  (l/3+3(c)*t)  tends,  as  e  -*  0+, 

to  a  cdf  4>^(t)  5  1—  4>^(t)  with  hazard  rate  x^(t)  (see  for  example 
[1]  for  discussions  of  the  concept  of  hazard  rate  in  its  usual  context 
of  reliability  theory),  which  we  call  I's  asymptotic  distribution  of 
time  of  (first)  fire.  Analogously,  x^Ct)  is  the  hazard  rate  of  ll's 
asymptotic  distribution  ^2(t)  of  time  of  first  fire. 

Hence,  if  T.  _  is  Player  I's  time  of  (first)  fire  under  mesh 

-  C’1  •' 

c,  (Tj  f.-l/3)/8(c)  is  seen  to  converge  i •'  law  to  the  cdf  with  hazard 
rate  x^,  and  an  entirely  analogous  argument  and  conclusion  hold  for 
Player  II.  Normalization  is  by  the  location  parameter  1/3  and  the  scale 
parameter  8(c),  the  latter  in  effect  stretching  the  shrinking  interval 
[p(c),  1/3]  into  the  unit  interval  J-1,0].  Of  course,  in  view  of 
Theorem  2  (page  42)  of  [7],  essentially  no  other  normalization  will 

A. 

work. 

Having  obtained  these  results,  we  can  now  look  back  on  their 
derivations  and  sec  why  Hypotheses  A  and  B  constitute  an  efficient  set 
of  assumptions.  We  do  this  by  attempting  to  alter  Hypotheses  A  and  B 
to  the  natural  alternative  Hypotheses  A’,  B’  and  1)  below,  equally 
plausible  in  the  light  of  the  computed  runs,  and  showing  that  we  are 
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driven  back  to  A  and  B. 

To  begin  with,  we  now  postulate,  with  0(c)  =  e/3(c)  and 
3(e)  -  1/3  -  p(c)., 

HYPOTHESIS  3(c)  -*■  o"*"  as  e  -*■  O"*",  which  implies  that 

P (t , c)  -  1/3, 

HYPOTHESIS  B*.  There  exists  a  function  v(t)  on  [-1,+  “')  such 

•J 

that  (u(t,c)  +  1/3) /6(c)  has  limit  -  v(t)  as  c  •+  0+, 

l1- 

and 

HYPOTHESIS  n.  There  exists  a  function  x^(t)  on  [-1,+  °°)  such 

that  x^(t ,c)/0(c)  has  limit  x^(t)  as  e  -*■  0+. 

Then  (10)  implies  the  existence  of  the  limit,  as  c  -*■  0+,  of 
(p(t,c)  -  l/3)/0(c)  on  [-1,+“)  ,  and  the  definition  of  p(t,c)  then 
implies^the  existence  of  the  limit,  say  C,  of  B(c)//e  as  e  0+, 

A 

which  in  turn  implies,  by  relation  (12),  the  existence  of  the  limit 
x2(t)  of  x2(t,e)/0(e)  on  [-1,+  °°)  as  C  -*•  0+. 

Now,  if  it  were  true  that  0=0,  then  (11)  and  (13)  would  say 
Xj(t)  “  9v(t)  =  x2(t)  for  all  t  >_  -1.  Furthermore,  we  could  arrive 
at  the  key  relation  (8)  with  the  left  member  replaced  by 

-  (v(t+0(c))  -  v(t) J/0(e)  +  lg(t+0(c),c)  -  g(t,c)]/02(c), 

and  there  seems  to  be  no  reasonable  way  to  utilize  this  relation  with¬ 
out  strengthening  Hypothesis  B’  to  Hypothesis  B,  in  which  case  the 

2 

differential  equation  (14)  becomes  v'(t)  =  9v  (t),  while  the  initial 
condition  (17)  becomes  v(-l)  =  0.  But  the  only  solution  of  this 
initial  value  problem  is'  v(t)  =  0.  'fhat  in  turn  forces  x^(t)  H  0, 
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corresponding  to  no  non-degenerate  cdf's  <5>^(t).  This  we  are  forced  to 
conclude  that  C  >  0,  which  amounts  to  strengthening  Hypothesis  A'  to 
Hypothesis  A. 

5.  Comparison  of  asymptotic  distributions.  Some  relationships 
between  the  cdf's  ♦  (t)  can  be  derived  under  hypotheses  weaker  than 
A  and  B.  If  we  assume,  again  based  on  the  pre-asymptotic  computer  runs, 

that  A  holds  but  that  B  is  weakened  to: 

*■  '  * 

B':  For  t  >  -1,  lirn  U-(t =  -  v(t)  exists, 

£~»0+ 

then  our  derivations  in  Section  3  obtain  the  existence  and  non-negativity 

of  the  limits  x  (t)  =  lim  x  (t,£)/0(e),  and  also  formulas  (11)  and  (13) 

e-0+ 

for.  Xj(t).  From  these  formulas  we  can  see  that  the  hazard  rates  are 
related  by'  -  .  x 

(29)  x2(t)  -  x  (t)  +  27C?t  . 

In  the  derivation  of  the  formulas 

(30)  ^(t)  *  exp(-  /^ixi(a)di) 
we  needed  A  and 

B":  Hypothesis  B'  holds  and  v(t)  is  continuous. 

Thus,  under  A  and  B",  relations  (29)  and  (30)  give 

(31)  ’  ^2(t)  ■=  expld2  (1-t2)]  *.(t), 

I . 

where  <S  ■  C/27/2.  This  implies  the  comparisons 
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(32)  <J^(t)  >  ^(t)  for  -  1  <  t  <  1  and  #  (t)  <  4>2(t)  for  t  >  1, 

which  is  to  say  that  $2  *s  more  highly  concentrated  about  the  normalized 
firing  time  l  than  is  •J'j* 

Hypotheses  A  and  B"  also  imply  that  the  densities 

(33)  fi(t)  3  -  4>^( t)  »  xi(t)  exp(-  /^xi(T)dT) 

»» 

are  defined  and  rton-negative.  However,  it  is  not  apparent  that  rigorous 
analysis  of  these  densities  can  be  done  without  invoking  the  full  force 
of  A  and  B.  Since  these  hypotheses  eventually  lead  to  formula  (21)  for 
v(t),  it  is  evident  that  the  analysis  of  the  densities  depends  critically 
on  the  properties  of  the  solutions  iHs)  of  problem  (22).  For  example, 
formulas  (11)  and  (21)  imply  that  x^t)  =  -  ln|i|i(s)(,  hence  that 

(34)  ^(t)  *  <Ks)/tM-6)  • 

i. 

As  is  shown  in  {9],  <Ks)  >  0  for  all  s  and  iKs)  0  as 

00 

s  ■  fit  •+  +°°,  so  that  /_j^(t)dt  *  +  °°,  and,  in  view  of  (33), 

(35)  /”1f1(t)dt  -  1; 

this  is  to  say  that,  even  though  I's  asymptotic  firing  time  distribution 
is  heavy-tailed  compared  to  that  of  II,  the  convergence  of  I's  normalized 

firing  time  distribution  is  complete. 

*'  00 

Relation  (29)  implies  that  I  ^x2(t)dt  »  +~  also.  Hence  (33) 
yields  as  well  that  /^^(Odt  *>  1. 

Inspection  of  figure  1,  based  on1, the  computer  runs,  suggests  the 
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results  that 

(36)  (i)  mode  (<J>^)  =  -1  and  (ii)  mode  (G^)  >  0. 

If  both  A  and  B  hold,  then  the  results  of  [9]  can  be  used  to  establish 

(i),  based  on  (34)  and  the  facts  that  iKs)  is  positive,  decreasing, 

and  convex,  and  also  to  derive  (ii)  from  a  special  relationship  between 

iHs)  and  the  Weber  parabolic  cylinder  functions  D  and  °2/3*  In 

**•  ’  ’ 
particular,  the  differential  equations  analysis  leads  to  an  approximation 

mode  (<t.j)  =  1.5058,  furnished  by  R.  J.  Lambert,  which,  as  figure' 1  suggests, 

is  in  accordance  with  the  computer  runs. 

No  amount  of  inspection,  however,  will  suggest  that 

(37)  (i)  mean  (G^)  <  00  and  (ii)  mean  (G^)  *  <*>  . 

A 

Under  A  and  B,  we  can  use  (31),  (34),  and  the  definition  of  G^(t)  to 
get 

(38)  *2(t)  -  exp(62(l-tz)]  iKs)/iK-5). 

The  facts  that  <l>  is  positive  and  decreasing  force  G^(t)  to  be  trapped 
between  0  and  exp[62 (1-t2) ] ,  whicluleads  to  (37(1)) .  The  infinite 
expectation  for  player  I  is  established  in  [9],  based  on  the  fact  that 
G^(t)  *  G(<5t)/iH-5)  behaves  like  t  as  t  -*■  «°  . 

6.  The  m  vs.  n  case.  We  conclude  with  a  brief  outline  of  how 
results  analogous  to  the  above  are  obtained  for  the  m  vs.  n  equal- 
accuracy  noisy  duel,  where  m  <  n.  The  details  of  the  derivations 


can  be  found  in  19]  and  (10]. 

The  analysis  begins  with  the  replacement  of  the  matrix  (2b)  by 
F  7 

F  (1-p)Vl,n-l[p+el  P  +  <1“P)Vl,n(«,+£j 

F  -  p  +  (l-p)um>n_1(p+c]  Um,nlp+C] 

M 

where  u.  (p+c]  is  the  value  of  the  i  vs.  j  game  at  the  grid  point, 

i.j  *• 

p+e.  In  particular,  this  value  is  0  if  i  =  j  and  is  -1  if  i  =  0. 

Hypotheses  A  and  B  are  modified  suitably,  again  based  on  computer 
runs  for  some  small  values  of  m  and  n.  Based  on  those  hypotheses,  we 
can  adapt  the  methods  of  Sections  3  and  4  above  to  derive  the  following 
formulas  for  the  asymptotic  distributions  ^  and  ^  of  the  players' 
normalize'd- times  of  first  fire. 

The  cdf  (t)  »  \Ks)/iHs  ),  where  s  -  5(t+n)  with  6  >  0,  -  n 
1  o 

is  the  value  of  t  >  -1  where  the  two  players’  hazard  rates  are  equal 
(which  is  zero  in  the  case  1  vs.  2),  s^  =  <5(-l+n),  and  i|>  is  the  solution 
of  the  Hormite  equation  _ 

i|i"(s)  -  2sip’(s)  -  4ai[(s)  ■  0,  a  =  m/2(nrt-n), 

% 

that  satisfies  the  initial  conditions 

«H0)  -  1,  <|»’(0)  -  -  2r(a+‘i)/r(a). 

The  cdf  ( t )  -  G(t)  ‘I’j(t),  where  G(t)  -  exp(s2  -  s2). 

Based  on  these  formulas  and  the  properties  of  the  function  ip 


derived  in  [9],  we  find  that 


a 


4>^('t)  >  ^(O  for  -1  <  t  <  l-2ri,  4>^(t)  <  ^(t)  for  c  >  1-2n 

lim  4>  (t)  *=  1 

t-*oo  1 

and  that  mode  (4>^)  =  -1  <  -  n  <  mode  (<£>2)  while  -1  <  mean  (4>2)  <  mean  (4^) , 
with  mean  =  +*>.  Thus,  we  have  the  inequalities  (la)  and  (lc) 

for  the  m  vs.  n  case  while  maintaining  the  weaker  player's  complete 
convergence  of  normalized  firing  time,  and  infinite  normalized  firing 
time  expectation. 
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